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Abstract 

Given a real objective function denned over the sym- 
metric group, a direct local-search algorithm is pro- 
posed, and its complexity is estimated. In particular 
for an n-dimensional unit vector we are interested in 
the permutation isometry that acts on this vector by 
mapping it into a cone of a given angle. 

1 Introduction 

In several areas of Game Theory and Oil Engineer- 
ing there appears an optimization problem defined 
over the group of permutations on an index set. The 
declination curves [1] on wells, for instance, are im- 
portant in order to develop exploitation strategies, 
among them Steam Assisted Gravitational Drainage 
(SAGD) processes [3J. The well known theories de- 
veloped by Von Neumann on Game Theory pose also 
optimization problems whose domain is the symmet- 
ric group. A rather complete survey of this kind of 
optimization procedures in Social Sciences is given 
in 0. 

A set of independent measurements of a certain 
parameter is realized as a vector in M" and when 
it is normalized (in the statistical sense: it is trans- 
lated by its average and elongated by the reciprocal 
of its variance) an unit vector with zero mean and 
unit variance is obtained. Then the problem consists 
in finding a permutation of its entries making a given 



angle with the original normalized vector. 

Here we sketch a local-search procedure, with cu- 
bic complexity with respect to the dimension of the 
original vector. 

2 Uniform distribution 

Let e, = (S{j)" =1 be the i-th vector in the canonical 
basis of R n , where 8ij is Kroenecker's delta. For any 
non-zero vector x = (x%, . . . ,x n ) = Yjj=i x j e j' ^ s 
mean is E(x) = ^ X)j=i x ji ^ s variance is Var(x) = 
~ yVj—l (xj — -E(x)) 2 and its normalization is 

zx= . 1 = = (x-E(x)l") (1) 
V Var(x) 

where 1™ is the vector in K™ with constant compo- 
nents 1. Clearly, E(z x ) — 0, Var(z x ) = 1 and z x is 
an unit vector, with respect to the Euclidean norm. 
In particular, for the identity vector 



(1,2,. 



its mean and variance are 



i=i 



(2) 
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and its normalization is z = y^_i z 7 -e,-, where, for ,„ . . ,.. ,„ ... , 

eachi<n: ° - (n + 1)) (2*0) - (» + 1)) 



t , Let us write v n = _ Prom eq. (TIT)]) it can 

be checked that the image of <!> is included in the 



j(n - lWn + 1) , ., discrete set 



(6) «n - < — 



Then 



fceZ, -v n < k < v n \ c [-1,1] (11) 



_ , (n - l)n(n + 1) Namely, $(5„) C P„, although the inclusion may 

' " s 3 be proper, for some values of n. For instance, by 

eq. (Uni), Vfc G Z, 

— G $(Sn) ^=> 

Let «S n be the permutation group of the set of indexes 

[l,n]. Then cardS^ = n\. A permutation matrix J\ 

P G R" xn has the form P = (d" ir(i) ). . for some 3?r G : fc = ^ (2j - (n + 1)) (2?r(j) - (n -+(1^) 



3 Optimization problem 



' i,j<n 

permutation n E S n . In this case, let us write P — 



3 = 1 



Ptt. This matrix determines a linear isometry P n : w , . , , 

17 Whenever n is an even integer, each summand 

R" — * R": for each y e M™, P^y) coincides with y , . / , , . , , / , 

, , ./ , ' |fp r\ ( v) ( )) odd (as the prod- 

up to the permutation tt, hence ||n(y)|| 2 = ||y|| 2 . uct of tw0 odd in t e gers), thus their addition will be 

For any non-zero vector x = £,=i Xj e, G R", eyen _ Hencej for fc oddj wg haye _k_ ^ 

by Schwartz's inequality, the inner product between „ , ,, . " x /cl > . 

, t , , • i. Besides, we can see that the image $ i„ is sym- 
the normalized vector and its corresponding image 



,, „, metric with respect to 0, namely, W G S n 

satisfies |(z x ,P„(z x )| < ||z x || 2 ||P w (z x )|| 2 = 1. Thus, 

for all 7r G S n , n 

fc = ^(2i-(n + l))(27r(j)-(n + l)) =^ 

-i<(z x ,p.(z x ))<i (8) fri 

n 

Let us consider the map -fc = £ (2j - (n + 1)) (27r„, o 7r(j) - (n + 1)) . (13) 

$x : S„ -> [-1, 1] , tt (z x , P 7r (z x )) (9) J=1 



For the identity permutation, vn = (1 2 ■ • • n) we Let us consider the following 
have $ x (7Ti) = (z x ,z x ) = 1 and for the permutation 
7To G SVj such that P 7ro (z x ) = z x we have $ x (7To) = 
-1. 



Problem 3.1 Given a non-zero x G R" and r G 

[-1,1] y>nd 



In particular, if x is the identity vector, n above is = _ 

the reverse permutation n n i = [n n — 1 ■■■ 1), namely " 



Vj, (27r n ,(j) - (n + 1)) = - (2j - (n+ 1)). Besides, 



From remark (|13[) . it is enough to assume r G [0, 1] 



from eq's (El) and ©, , . ., ' ., , 

i i— v v—in when x is the identity vector 

3 

(n- l)n(n+ 1) 



3 Problem 13.11 can be solved by a direct local-search 

z i-f7r(z)) — , — 7\TZ7 T 7\ ' algorithm. Let us consider the elements of S n as 



2 



nodes and let us consider edges of two kinds: 

(■K,p)eE T 3i,j G {l,...,n} : p = 

(ir,p)eE c 3i G {1, . . . ,n} : p={i,i')w 

with i = i + 1 if i < n & 71 = 1 

Let S nT = (S n ,E T ) and S n c = (SW, E c ). The graph 
S„t contains as edges all pairs of permutations that 
differ by an arbitrary transposition, and the graph 
S„c contains as edges all pairs of permutations that 
differ by a transposition of two consecutive indexes. 
S n T is a regular graph of n\ nodes of degree (™) = 
\n[n — 1). If we consider any edge as an edge of unit 
length, then the diameter of S n T is d n T = n — 1 = 
0(n). Instead, S n c is a regular graph of n\ nodes of 
degree n and its diameter is 

dnC — -r(n ~~ (timod2))(n + (nmod2)) = 0(n 2 ). 



Input. The dimension n 6 N and the objective map 
F:S„->R + . 

Output. The permutation it G S n such that ir — 
Arg Mm peSn F(p). 

1. Initialize n c := iri (the identity permutation); v c := 

2. repeat 

(a) ir d := 7r c ; v d := v c ; 

(Jo) let 7r c be the i?-neighbor of TTd minimizing the 
objective map F in the ^-neighborhood of TTd] 

(c) v c := F(tt c ); 

(d) if v c > v d then ir c := ir d ; 

until 7r c == ir d ; 

3. Output 7T := 7T C . 



4 Solution procedure 

Let F : 5„ — * M + be a positive-real valued map de- 
fined over the symmetric group. Let S„ = (S n ,E) 
be any of S„t or S„c- A local-search algorithm to 
minimize F is sketched at table Q] 

In order to solve problem 13.11 for any given x G 
K™ — {0}, r G [0,1], the objective map F xr : tt > 
|r — $ x (7r)| is considered, where $ x is given by eq. 
which in turn involves the normalization z x of vector 
x, given by eq. ([5]) for the identity vector. 

If S„ = S n T, then the graph is regular of degree 
n(n — l)/2 and diameter n — 1. Thus step 2.(b) at 
pseudocode in table Q] entails n(n — l)/2 evaluations 
of the objective map, and whole cycle 2. is repeated 
at most Ti—l times. Thus, in this case the worst case 
entails n(n — l) 2 /2 evaluations of the objective map. 

Instead, if S n = S n c, then the graph is regular 
of degree n and diameter d n c- Thus step 2.(b) at 
pseudocode in table Q] entails n evaluations of the 
objective map, and whole cycle 2. is repeated at most 
d n c times. Thus, in this case the worst case entails 
ndnC evaluations of the objective map. Since 



Table 1: Local-search algorithm on the symmetric 
group. 

we see that using S n T is twice as expensive than using 
SnC, and this last entails 0(n 3 ) evaluations of the 
objective map. 
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